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Abstract 

We study the construction of the collective-coordinate manifold in the baryon 
number two sector of the Skyrme model. To that end we use techniques of adia- 
batic large amplitude collective motion, which treat potential and kinetic energy 
on an equal footing. In this paper the starting point is the Ansatz proposed by 
Atiyah and Manton (Phys. Lett. 438B, 222 (1989)), which allows a study of the 
dynamics using a finite and small number of variables. From these variables we 
choose a subset of collective ones. We then study the behavior of inertial pa- 
rameters along parts of the collective manifold, and study the dynamical parts 
of the interaction. 

1 Introduction 

In this paper we study an attempt to model the nuclear force from an underlying 
effective field theory, that mimics some aspects of the "correct" theory of strong inter- 
actions, namely QCD. This is in sharp contrast to the standard effective models of the 
nuclear interaction, which have been derived mainly through a description in terms 
of (fictitious) meson exchanges. In general they have received remarkably little direct 
input from QCD (see Ref. U for a discussion of several approaches). Nevertheless such 
an approach is both a quantitative and qualitative success. 

Another possible approach, which has only been pursued in recent years, is the use 
of models inspired by the large N c limit of QCD pi |3|. The prime example is the 
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Skyrme model |4[], a non-linear theory of interacting pions, where baryons appear as 
topological defects - and this implements confinement. The model was conceived and 
forgotten in the early sixties, and was revived only in the early eighties by Witten 
and collaborators ||. It has since been used successfully to describe the properties 
of baryons. Of course it was soon realized that one could use the model to describe 
baryon-baryon interactions as well || [7]] . Actually even this idea seems to have been 
discovered by Skyrme ||. 

The baryon-baryon interaction derived in this model was plagued by a lack of 
central attraction. In the early works the interaction was derived using the product 
approximation, which assumes that the two Skyrmions do not deform as they approach. 
Furthermore one ignored the mixing between the nucleon and its excited states, which 
gives an attractive contribution to the central force. This question has recently been 
investigated in some detail, and it was found that the reintroduction of a finite number 
of colors (A c ), together with the use of a Born-Oppenheimer approximation is necessary 
to describe this part correctly |9|, [UJ . 

Apart from the static interaction (i.e., the part of the energy independent of ve- 
locities) there is a dynamical structure of the Hamiltonian as described by the kinetic 
energy. The inertial parameters are in general position dependent. If we use the fact 
that the mass tensor can be used as a natural metric on the configuration manifold, 
this shows that that coordinate space is curved. Actually little is know about these 
parameters. The only calculations have been performed using the product ansatz, and 
are thus only good for large separations. One of the important parts of the NN force, 
the spin-orbit interaction, derives purely from these kinetic terms. One finds that, even 
when the separation is large, and the product Ansatz is supposed to work, the sign 
of the spin-orbit force 0, O, O, IT3, O] does not agree with phenomenology. This 



problem has not yet been resolved. 

The potential and kinetic energies can be derived from a knowledge of the solu- 
tions of the Skyrme model for finite separations. Unfortunately no such solutions exist. 
The best calculations up to date |1| calculate a set of fields obtained by imposing a 
constraint that is chosen in advance, and not determined self-consistently. As an alter- 



native there exist the set of parametrized fields introduced by Atiyah and Manton [17 



These authors use an interesting technique to obtain a field of given topological quan- 
tum number from a parametrized instanton field. Its advantage is that the instantons, 
which are specified by a finite number of parameters, induce an effective Lagrangian 
for the Skyrme model depending on the instanton parameters. One can attempt to 
choose combinations among those that describe the collective manifold. As a first step 
towards this goal a calculation, similar in nature to that of Walhout and Wambach, of 



the static potential was performed by Hosaka et al ||18|| . This gave reasonable answers. 

Since one can perform exact numerical solutions for the same problem, the greater 
promise of this Ansatz lies in the understanding of the geometrical structure of the 
collective manifold, which is closely tied to the kinetic energy. Atiyah and Manton 
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have recently [|19[] analyzed some aspects of this problem analytically Several of the 
more quantitative aspects of this approach require numerical calculations, and cannot 
be resolved by analytical means. This paper is an attempt to study such an approach. 
There is a related idea to use the Atiyah-Manton Ansatz to study peripheral scattering 
2"Ufl. Our work has no overlap with this approach, since we shall concentrate on short 



and intermediate distances. We believe that the long distance behavior can be best 
described by the product approximation, maybe suitably symmetrized. 

The determination of collective coordinates in classical mechanics, or - equivalently 
- time- dependent mean field theories, has been studied extensively over the past 20 
years (see Ref. |21[ for a discussion). As discussed in Ref. |H] the theory of adiabatic 
large amplitude collective motion (ALACM) is developed well enough to construct 
algorithms that are singularly well suited for this task. In this note we shall apply the 
AM Ansatz, and mix it with the notions of ALACM in order to study the structure of 
the adiabatic manifold. Similar calculations can be - and should be - performed by 
putting the Skyrme model on a grid. Since such a task is much more complex, both 
numerically and analytically, it is much harder to gain insight in that way. 

The paper is organized as follows. In Sec. Q we succinctly introduce the Skyrme 
model as well as the product Ansatz and the Atiyah Manton Ansatz. In Sec. |3] we 
give a short introduction to the theory of Large Amplitude Collective Motion, as far as 
relevant for the current paper. In Sec. f| we introduce the numerical techniques used 
in our calculations. In Sec. |5| we discuss the results of these calculations. Finally, in 
Sec. H we gi ve a discussion and outlook. In three appendices we discuss which classes 
of solutions have a reflection symmetry, we give a detailed analysis of the fluctuations 
around the B = 2 hedgehog and we discuss the algorithm used in obtaining self- 
consistent solutions for the LACM equations. 



2 Skyrme model and Atiyah-Manton Ansatz 
2.1 The model 

The "standard" Skyrme model is based on the non-linear sigma model, extended by 
a quartic interaction term and a pion mass term. The Atiyah-Manton Ansatz can 
only be used in the (chiral) limit of zero pion mass, where the model is defined by the 
Lagrange density 

C = f l T r{d^Jd^) + -L^Tr^d^U, d v U][tf d»U, tfd u U], (1) 
4 dZg p 

where U is a unitary two-by-two matrix-valued field satisfying the boundary condition 
U = 1 at infinity. As has been discussed many times before, this model has a topolog- 
ically conserved quantum current. The charge of this current is identified with baryon 
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number B. Finally, one can also separate the U field in a a and pion field, 



U = ^-(<r + in-T). (2) 

JTT 

On occasion it may be useful to look at the vector pion field. 

If we rescale the units of time and length, x — > x/(gf n ), (the so-called Skyrme 
units), the Lagrange density takes on the slightly more convenient form 

C = O-Trid^Ud^) + —Tr[U^d^U, d„U][tf d^U, U ] d v U]\ , (3) 
9p 32 / 

where f n /g p is the Skyrme unit of energy^ Finally the Skyrme Lagrangian can easily 
be reformulated in terms of the Sugawara variables (Lie-algebra valued currents) L, 



L, = U- 1 d,U = il«T a , (4) 



and we have 



r = -7 a p + - 



1 r 



{iT a f - i a JT h i vh ] (5) 



2.2 Structure of the B = 2 manifold 

Preliminary investigations of the structure of the B = 2 manifold have been performed 
in great detail using the product Ansatz, and are discussed in the reviews p2| , p3| . 

Let us recapitulate the necessary details. For baryon number one the basic solution 
is the hedgehog, 

U = exp(zr • rf(r)), (6) 

with /(0) = 7r and /(oo) = 0. The baryon density has spherical symmetry, and 
the pion field points radially outward at each point of this surface. Actually one can 
perform a constant isorotation on the quantity U, AUA\ without changing the energy. 
This changes the direction of the pion field, and thus "grooms" the hedgehog. 

One can use these solutions to construct the so-called product Ansatz. This uses 
the fact that the product of two U fields each with baryon number one, with arbitrary 
centers and grooming, has baryon number two. As Skyrme already understood, this 
Ansatz is only good at large separations. At smaller separations the two U fields 
no longer commute, and an asymmetry between the treatment of the two hedgehogs 
appears. Furthermore, the bound state in the B = 2 sector of the model, named 
"donut" after the toroidal symmetry of the baryon density |24|, [27f , can not be 

described by the product Ansatz. 



^^There are many inequivalent definitions of the Skyrme units in the literature, depending on the 
definition of f w . In our definition the S.U. of energy is equivalent to 11.18 MeV (for the "standard" 
values of /„ = 54.1 MeV and g = 4.84), and the S.U. of length is 0.754 fm. 
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attractive 




hedgehog 

Figure 1: A schematic representation of the three channels, as they would appear in the 
limit of large separations. The spheres represent a surface of constant baryon density, 
and the arrows show the pion field on this surface. Here we have identified space and 
isotopic space. 



One way out is to solve the Skyrme model numerically on a grid, with imposition 
of the separation of the two Skyrmions. This approach has been taken by Wambach 



and his collaborators ||T6| , |28[| . Actually they never solve for all possible states of two 
Skyrmions, but limit themselves to a subset of configurations with definite symmetries 
under reflections, called the attractive, repulsive and hedgehog channels. These con- 
figurations are represented schematically, using the product Ansatz solutions in which 
they go over for large separations, in Fig. 11} (A similar study using the Atiyah-Manton 
Ansatz was performed in Ref. |T8| .) 

The three channels are first of all the hedgehog channel, where we perform no 
rotation on the two hedgehogs. All states have the same reflection symmetries as 
a single baryon-number two hedgehog. Second is the repulsive channel, where we 
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perform a 90° rotation of each hedgehog (180° relative angle) around the axis through 
the centers of the hedgehog. This is called repulsive since there is a strong repulsion 
in this channel. Finally there is the attractive channel where each hedgehog is rotated 
by 90° degrees about an axis orthogonal to the line connecting the two centers. 

From the caricatures of the pion field drawn here one can easily see that the fields 
have a certain symmetry under reflection in each of the three coordinate planes. These 
symmetries are discussed in greater detail in Appendix [A[ The idea is that these 
symmetries endow a special importance to these configurations, since they are at least 
stationary against small fluctuations breaking the symmetries. 



2.3 Atiyah-Manton Ansatz 



As discussed in |19| one can derive a Skyrme field from an instanton field by integrating 
the time component of the gauge potential, 

U(£) = Cs(Pexp J°° -A A {x,t)db\cK (7) 

Here S is a constant matrix, chosen such that U decays to 1 at infinity, and C describes 
an overall grooming. For the current work, where we shall only consider C near the 
identity, it is convenient to parametrize 

C = exp(ff -0). (8) 

For the Jackiw-Nohl-Rebbi (JNR) instanton of charge k = 2 we have pSfl 

A 4 (x,t) = - r, P = 2^1 V~V 9 

2 p ^ \x - Xi\ 

and we should use S = —I, to obtain a field of baryon number B = 2. 

In order to solve for the AM value of U we convert the integral (0) to the solution 
of a differential equation. First introduce 

U(x,t) = CS^PexpJ' -A A (x,t)dt^C j . (10) 

This function satisfies the differential equation 

d T U(x, r) = -U(x, t)A±(x, r), (11) 

with initial condition U(x, — oo) = S. The function U(x) is obtained as the limit for 
r — > oo of U. Actually we prefer to work with the Sugawara variables L M , Eq. (|j), 
so that the Lagrange density does not contain explicit derivatives. These can also be 
calculated as the large r limit of a quantity L M , which satisfies the differential equation 

d T L^(x, r) = [A 4 (x, r), L^x, r)] - d^A^f , r). (12) 
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Here the boundary condition is L^(x, — oo) = 0. By differentiating the differential 
equations (12 ) we can obtain expressions for derivatives of with respect to the 
instanton parameters A^ and Xi, which will be needed later. 

The field U defined in Eq. ([?]) does not have an explicit time dependence. There is 
an implicit dependence, due to a possible variation of the instanton parameters as well 
as the unitary matrix C with time. Let us denote the parameters {A;, Xi, 9} collectively 
by £. We then have 

L = i a U ] d ia U = i a L, a . (13) 
If we substitute this in the Lagrangian we obtain the form 



T — V, 



with 



and 



T 



B, 



a/3 



E M + 2 



E^ a ) -e(eW) 

i i,a I ij \ a / 



E W E #i -EE W% E #1 



(14) 



(15) 



(16) 
(17) 



The Lagrangian is quadratic in the time-derivatives due to the special nature of 
the Skyrme model. This also means that we have broken the Lorentz invariance of the 
original equations, so that (p~3|) can only be used to describe adiabatic (small velocity) 
motion. 

In the Atiyah-Manton Ansatz we have thus replaced the general matrix U(x), 
with an infinite number of parameters, by a form parametrized by 18 parameters, 
CU(x\£)C*. (At each point we can extract at most 16 dynamical parameters, we can 
always form at least 2 spurious combinations. For the channels with additional sym- 
metry discussed in this paper, we always find one more spurious combination.) Even 
these 15 or 16 dynamical variables are too much to describe the collective motion. 
We wish to study the behavior of the slowest modes among those appearing in (|T3"D. 
The techniques for such an approach, in the case that the Hamiltonian is quadratic in 



momenta, are well-developed plf . In the next section we will recapitulate one of the 
possible mechanisms of solution, that will prove convenient in the current context. 
Let us now enumerate those configurations, that will be relevant for the rest of the 



discussion. Here we closely follow [Fj], 
1. Product Ansatz. 

The choice A2 Ai, A3 describes a product Ansatz solution with centers at X\ 
and A 3 . In this limit one finds 



U(x) ss GlU^x - X 1 )G 1 G\U 1 {x - X 3 )G. 
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Figure 2: A schematic representations of the pole positions that describe the three 
channels in the Atiyah Manton Ansatz. The pole configuration a) represents a B = 2 
hedgehog, b) states in the hedgehog-hedgehog channel. The triangle in c) represents 
the attractive states, and d) the repulsive states. 



where the iso-rotational matrices G are given by 



(Ti - T 2 ) + i(Xi - A 2 ) • f 



2. B = 2 hedgehog. 

The hedgehog solution is found for the choice Ai = A 3 = 1, A 2 = 32, Xi = 0, 
Ti = -T 3 = 9.5 and T 2 = 0. 

3. The hedgehog channel. 

In general states in the hedgehog channel are described by a pole configuration 
with all three poles on the x-axis, Ai = A3 = 1, A 2 = A, An = — A31 = —A, 
A 2 i = 0, Aj( 23 ) = 0, Ti = 0. The problem is that, even though the states 
in the hedgehog channel do have the same symmetry as the hedgehog, there 
no continuous path from one to the other conserving the symmetry. The only 
common point, where we let the three poles coincide in the origin, corresponds 
to a state with B = 0. We shall turn to this problem in a future section, when 
we study the hedgehog numerically. 

4. The attractive channel. 

Here we have Ai = A 3 = 1, A 2 = A, A n = -A 31 = -A, A 12 = A 32 = Y, A 22 free, 
all remaining components 0. Actually it was argued in [I8| that a simple relation 
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between the parameter A2 and the shape of the triangle (see below) describes the 
attractive channel best. 



5. The repulsive channel. 

Here we have the three poles aligned along the x-axis, Ai = A 3 = 1, A 2 = A, 
Xu = —X31 = —X, T\ = T 3 = T, T 2 free, all remaining components 0. Again 
the authors of Ref. [|ll| used the same relation as in the attractive channel between 
the shape of the triangle and A. 



Let us use the ideas of [18j to study the configurations in the three channels. The 
authors of this paper use what can be called a turning-point approximation: At each 
point of the path they only consider the potential energy as if it were a stationary 
quantity. Of course this is only true if we tune the total energy at which we work to 
be equal to the potential energy. We thus end up in a turning point, since the kinetic 
energy has to be zero! Furthermore we can then use a virial argument to show that 
the quadratic and quartic terms in the Skyrme Lagrangian contribute in equal parts 
to the energy. 

The formalism in the next section tries to improve on such an approach by including 
information about the kinetic energy. As as a first guess the turning-point approach 
is fast and cheap. Let us look at the attractive and repulsive channel, and study the 
dependence of the energy on the form and weights of the triangles. The only remaining 
parameter, the size of the triangle, can be determined from a virial argument, since 
E 2 = E4 at a stationary point. The hedgehog-hedgehog channel is of no special 
interest since it is specified by only two parameters, one of which follows from the 
virial argument. The remaining parameter then specifies the separation of the two 
solitons. 

In order to study the energy we parametrize the triangles by introducing a second 



parameter At- In Ref. [18 it was assumed that this parameter equals \w, the weight 



A2. The parametrization takes the form 

Ai = A 3 = l, 
A2 = Avf, 



Xu = -D/y/(l - 1/(1 + A T ) 2 ), 

A31 = D/\/(l-l/(l + A T ) 2 ), 

X l2 = -D/(1 + X T ), 

X 32 = -D/(1+X T ), 

X 22 = D (20) 

for the attractive channel. For the repulsive case X 2 should be exchanged with T. 
We now study the energy as a function of At and \w for constant separation. The 
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Figure 3: The energy at constant separation in the attractive channel as a func- 
tion of log(Av^/Ax). From bottom to top the curves correspond to values of R of 
1.2,1.3,1.4,1.5,1.7,2.0,2.5,3.0 S.U.. Dashed lines have A r < 1 and solid lines A r > 1. 



definition of separation is somewhat arbitrary. We use a rotationally invariant defi- 
nition through the size of the quadrupole moment of the baryon-number distribution 

({f(x)) = fd*xB(x)f(x)) 



R 4 /4 



(x 2 ) 2 + (yY + (zY 

-(x 2 )(y 2 )-(x 2 )(z 2 ) 



(y 2 )(z 2 ) + 3((xy) 2 +(xz) 2 + (yz)' 



(21) 



The definition in Ref. |L8| is based on the 11 component of the quadrupole tensor. As 
all sensible definitions should, both reduce to the separation between the centers of 
two isolated Skyrmions for large R. 

In Fig. |3] we show the results for the attractive channel. Here the curves represent 
the energy at constant separation. It appears to be convenient to plot them as a 
function of log(Xw/^T)- We find two solutions of minimal energy for given R, both for 
Xw = At- This is due to the known inversion symmetry under A — > 1/A, where A > 1 
corresponds to separation along the x-axis, and A < 1 to separation along the |/-axis. 
(There may be a hint in our calculation that slightly unequal values of Xw and Xt may 
actually be the optimal solution, but that is probably due to a numerical inaccuracy.) 
A careful numerical study of the pion field suggests that the only solution with the 
requested mirror symmetries is the one with equal A's. Thus the valley in the energy 
for fixed separation corresponds to the equal values of A as used in Ref. ]T8 ]. Of much 
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Figure 4: The energy at constant separation in the repulsive channel as a function of 
\og(X w /X T ). The value of R is 1.3 S.U.. 

more interest is the repulsive channel, Fig. ^. Again a numerical study shows that only 
for \y/ = Xt we find the correct mirror symmetries for the repulsive channel. We have 
chosen to look at a small value of R, since only there the energy depends strongly on 
\w/Xt, even though the same trend appears to persist for all R. Here we find that the 
energy is actually a maximum in the configuration with equal A's (and thus having the 
symmetries). This is not unexpected, since one would like to see that in the repulsive 
channel only those solutions most repulsive in nature appear. 




3 Large amplitude collective motion 



The approach used in selecting collective coordinates is discussed in detail in Ref. |2T 
In this section we shall discuss the most important ingredients of that work, as well as 
the trivial transformation from a Hamiltonian to a Lagrangian formalism. 

3.1 General Principles 

Given a classical Lagrangian of a general quadratic form in time derivatives, 

L = \tB a ^)i P -V(0, (22) 

we look for those motions that are approximately decoupled. A decoupled surface £ is 
defined by the property that if we start the time-evolution from a point on the surface 
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(i.e., £ G £), with velocity £ tangential to this surface, the motion remains on the 
manifold for all times. Since a slight perturbation of the Hamiltonian will destroy such 
a manifold, one sees that only very special models will have exact decoupling. 

For that reason we look for approximate solutions. First of all, we select a subset of 
the exact decoupling conditions, that allows us to construct an algorithm to calculate 
a candidate manifold. Secondly we construct a decoupling measure, that allows us to 
gauge the quality of decoupling. This tells us whether it makes sense to consider the 
limited dynamics on the manifold. 

Let us start with the first problem, which is the more complicated of the two. As is 
argued in Ref. |?T), the mass matrix B a p actually plays the role of metric on the space 



q. We thus take over the standard notations of general relativity. In particular B with 
both indices raises is the inverse mass matrix. The geometric structure of the problem 
is especially clear in in one of the algorithms useful in tackling the decoupling prob- 
lem, the local harmonic approximation. This consists of the solution to the equations 
(the subscript ", a" is a shorthand for the partial derivative d/d£, a , "; a" denotes the 
covariant derivative, to be defined below) 



N 

V a = E AW /F, (23) 

i 

V^B^f® = (oo^) 2 f^. (24) 

The first of these equations states that the force is parallel to an approximate set of tan- 
gent vectors to E, while the second equation is a local harmonic approximation (LHA) 
to the motion, which defines the approximate tangent vectors. Of great importance is 
the use of the covariant second derivative of V, 

= \B lS {Bsp, a + Bg a ,p — B a f3 } s)- (25) 

We obtain the ordinary local harmonic approximation only at extrema of the potential 
energy (V a = 0). At these points Eq. (^) is automatically satisfied. This allows for a 
path tracking approach, where we bootstrap from a stationary point. 

Let us expand on this point, to clarify our approach. First note that since Eq. (|24"D 
is not a symmetric eigenvalue problem, we have independent left eigenvectors g". With 
the standard normalization condition 

9?Ji = & (26) 

these can be interpreted as defining the derivative of the old coordinates with respect 
to the new set which exhibits decoupling. The eigenvectors / then define the derivative 
of new with respect to the old coordinates. Suppose now that we start from a stable 
solution, in the case of the B = 2 Skyrme model that could be the donut or the 
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hedgehog. The direction in which we search are actually the least stable ones (the 
ones with the lowest harmonic frequency squared, (cu^) 2 ) in Eq. fl24[) . Suppose that 
we follow a single mode, as we will do in the next section. We then try to track a 
one-dimensional manifold by looking for a solution such that the RPA mode g\ is an 
approximate tangent vector to the path. We thus take a small step in this direction, 

r = £o+<i- (27) 

In general this point does not lie on the path, i.e., it does not exactly satisfy the force 
condition (|2"3"|), so that we have to look for a solution near this point. This is a non- 
linear problem, and can be solved in a variety of ways. In this work we use a very 
simple iterative method based on an approximate linearization of the force equations, as 
discussed in Appendix 0. Once we have located a second point on this path, we proceed 
with the same steps sketched above. For the problem of the Skyrmion, where we need 
to consider at least 6-dimensional dynamics to restore all the broken symmetries, we 
should actually calculate a many dimensional surface. We shall not pursue such an 
approach here, since it leads to significant topological and numerical complications, 
and shall concentrate on a small subset of the collective manifold instead. 

Once we have found the coordinate surface, we define coordinates functions Q 
decsribing it, and we introduce 

V{Q) = V(£(Q)), (28) 

and a collective mass B, equal to 

B^ = glB^g% (29) 

As usual for any linear eigenvalue problem, the scale of the eigenvectors is undetermined 
in Eq. (0). This freedom corresponds to a rescaling of the coordinates, which is always 
possible. Tensors are also transformed under such a rescaling. Of particular interest is 
the mass along the collective path, which can either be evaluated as 

B n = gftBfgfc, (30) 

or as 

Bn = -Qj-B^. (31) 

This last definition is based on the exact tangents to the decoupled manifold. Thus 
these definitions are only equal for exact decoupling. The normalized difference can 
be used to define a decoupling measure. If we normalize g such that Bu = 1, we find 



that this decoupling measure can be cast in the suggestive form p0 



D = YsW/dq 1 ) 2 . (32) 

i>2 
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Here one can approximate the change in the new coordinates, dq l , by a finite difference, 



dq l « f:JC, 



(33) 



which indeed defines a change in the new coordinates. It is obvious that if all coordi- 
nates are independent from the first coordinate we have exact decoupling. The steps 
dq 1 can then be used to define a collective coordinate, 



3.2 Redundant coordinates 

One important difference between the current problem and all the previous ones studied 
is the appearance of redundant coordinates. This is easily treated using the local har- 
monic approximation. The manifestation of redundant coordinates (that may change 
from point to point) is a zero eigenvalue of the mass matrix B, which is thus no longer 
invertible. But to calculate the covariant derivative we need this inverse! The solution 
is to only be perform the matrix inversion in the subspace with non-zero eigenvalues, 
the resulting T correctly eliminates these non-physical eigenvectors, including their 
derivatives, from the LHA equations ([M|). Here it is thus of utmost importance to 
use covariant derivatives, since the affine connection contains information about the 
change of the unphysical modes with a change of coordinates. For the present model, 
where the extra degrees of freedom can be eliminated explicitly, this procedure can 
easily be shown to give correct results. We have chosen to use the full set of degrees of 
freedom since it provides a sensitive test for the numerical accuracy and the correctness 
of our program. 

4 Numerical techniques 

From the previous sections one may notice that we need the spatial integration of the 
Lagrange density, and the second derivative of the potential energy as well as the first 
derivative of the mass. Since the potential energy depends on the current we need 
up to the second derivative of with respect to £ a . For the derivative of the mass, 
where the mass itself depends on If and l a a , we need additionally the derivative of l^ a . 
We have chosen to write exact differential equation for all these quantities, so as not to 
loose to much accuracy on the outset. Thus we need to solve, at each point in space, a 
set of 1269 coupled differential equations. This is solved by first going to the variable 
t = 7i / 2 arctan(t), and integrating from —1 to 1, using a highly efficient Runge-Kutta 
routine |H[ to solve the initial value problem. 

The spatial integral is decomposed into a radial integral and an integral over the 
surface of a sphere, which we deform to an ellipsoidal shape by multiplying each of 



Q = E dq 1 - 



(34) 



path 



14 



the major axes with a constant. The integral over the surface of the sphere was 
performed with an efficient non-product formula (see Ref. |32 for general references to 
such formulas), which also has at least tetrahedral symmetry. The relevant methods 
for the sphere are discussed in great detail in the Russian mathematical literature 
53L K3, RHj. For the work reported here we use both an 85 point 15th order and a 110 



point 17th order formula. These have proven adequate for the task. 

The remaining radial integral is transformed to the finite interval (—1, 1) by the 
mapping p = (r — a)/(r + a). We then perform a Gauss- Jacobi integration over p, 
using the weight function (1 — p) 2 (l + p) 2 . This takes into account the exact large r 
behavior of the slowest decaying integrand, which is the quadratic contribution to the 
potential energy, as well as the r 2 behavior near the origin. Judging from the values 
we found for the baryon number, this procedure is highly accurate. (We found a value 
of 2 in at least six significant digits.) Since the integrands are much more spread out 
the accuracy of the RPA matrices is somewhat less, but we feel that we can control 
the numerics. 



5 Results 

5.1 Attractive channel 

We start the application of the machinery discussed in the previous sections by study- 
ing the fluctuation around the minimum energy solution, the donut. In the AM Ansatz 
this state can be described by a configuration where all three poles have equal weight 
Aj = 1, and the poles themselves are located in an equilateral triangle, see Fig. ||c. 
We first study the harmonic modes around the donut. The eigenvalues and the eigen- 
functions (g a , the left eigenvectors of VB^ 1 (Eq. (0), that describe the change of 
the old coordinates as we follow each of the modes over an infinitesimal distance), are 
represented graphically in Fig. [5]. 

There are eight zero-modes, out of a maximum of nine. These correspond to the 
translations in space, rotations in space and rotations in iso-space. One disappears 
due to the special cylindrical symmetry, 2J 3 + J 3 = 0, of the donut. Let us discuss 
the zero-energy part of the spectrum first. The first and second mode form a pair, 
describing a rotation around the y and x axis, respectively. Since a rotation of the 
poles around these axes also generates an overall isorotation, these modes mix with 
an overall iso-rotation around the same axis to compensate for this effect. The third 
mode describes a translation in the ^-direction. The fourth and fifth modes describe 
iso-rotations around the y and x axis. The sixth mode describes the effect of the 
remaining (iso)rotational symmetry around the z-axis, ^3 — 2J3. Modes seven and 
eight are the translations along the x and y axes (slightly mixed). 

In table [I] we list the masses of the zero-modes, when normalized in such a way that 
they describe the effect of the standard generators. One should notice that the mass 
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Figure 5: A schematic representations of the eigenvectors of the RPA problem around 
the donut. The arrows in each box denote the direction of the given component of 
the eigenvector. The triangle on the left gives the x, y components for each of the 
three poles. The triangle itself denotes the position of the poles. The four rows of 
three points on the right denote, from bottom to top, the z components, the time 
components and the change in A, for (from left to right) poles 1,2,3. The upper line 
gives the change in 9i (Eq. (§)). 



Table 1: Masses of zero modes. 



type of mode 


mass (S.U.) 


translation 


X 


67.558 


translation 


y 


67.558 


translation 


z 


76.652 


rotation 


X 


86.964 


rotation 


y 


86.964 


rotation 


z 


130.342 


isorotation 


X 


57.258 


isorotation 


y 


57.258 


isorotation 


z 


32.585 


off-diagonal 


z 


65.171 
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Figure 6: The pion field in the plane obtained by following the ninth (lowest two 
pictures) mode and the tenth mode (upper two) by a very small distance. We follow 
the modes in both positive and negative directions. 



for translation is not completely isotropic. This is of course inconsistent with Galilei 
invariance, and appears to be incorrect. One can show that if the AM-donut were the 
minimum energy solution of the Skyrme equations this problem would not appear. It 
is due to the fact that the current Ansatz is not an exact stationary solution, and our 
covariant derivative, evaluated using the same Ansatz, can not correct for this effect. 
It shows some of the limitations of the present approach. 

Now let us look at the finite energy modes. The ninth and tenth modes are the 
lowest energy ones, and should be the modes we are looking for. The simplest way to 
understand the meaning of these modes is to look at a plot of the pion field in the xy 
plane. In this plane the pion field has no z component, so we can plot it as an arrow 
in 2D space. Such a plot can be found in Fig. where we have followed both these 
modes by adding and subtracting a small bit of the relevant mode vectors to the poles 
specifying the donut. The idea is that the pion field is only zero at the center of the 
donut, and that any perturbation will be clear from a change of the position of the 
zeroes. This is clearly the case! As one can see in the lower two figures, the ninth 
mode corresponds to the donut coming apart along either the x or y axis, the axis 
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depending on the sign of the mode eigenvector. The other mode corresponds to the 
system coming apart along an axis making an angle of 45° with the x and y axis. The 
reason is of course that the donut also has a reflection symmetry with respect to these 
axes! Actually the change in the pion field is trivial for these modes; it is a vector field 
parallel to the x or y axis, respectively. 

The eleventh and twelfth modes are somewhat harder to understand. They look 
somewhat similar to the previous two modes, if we rotate the coordinate system by 
45 degrees. They are very different, however! The pion field changes in a complex 
and much less symmetric way than in the case of the previous modes. In general the 
trend seems to suggest that these modes may be those that connect to the hedgehog- 
hedgehog channel. The thirteenth mode obviously describes the breathing mode of 
the donut. Finally the last two modes describe a rotations in the time plane. Naively, 
as discussed in Appendix |A] this is another class of modes that are compatible with 
the reflection symmetries. These modes are pushed to a high energy due to a strong 
mixing with the isorotation "C" degrees of freedom. One might be tempted - as the 
author was - to forget about these global isorotations. This is not even approximately 
right, however, since without C these highest two RPA modes become the lowest in 
energy! For large distances a similar rotation in the time plane on two widely separated 
hedgehogs corresponds to a zero mode, i.e., in that case this mode does not change the 
energy of the system. 

One of the main goals of the present work was to map out a path in configuration 
space where one follows the lowest mode, starting from the donut. This appears to be 
rather unproblematic (cf. discussions below), and there is no major obstacle. Some of 
the results are presented in Figs. [7|-[T0|. 

In Fig. [7] we exhibit the diagonal inertial parameters for rotations, isorotations 
and translations (i.e., the mass) as a function of the coordinate Q. These show the 
expected (and not very exciting) behavior. The mass is almost constant since the 
energy is almost constant, and the breaking of Galilei invariance remains a minor 
nuisance. The moment of inertia for rotations increases drastically as Q increases, for 
rotations orthogonal to the axis of separation (x-axis). Parallel to this axis it decreases 
as expected. There is only one off-diagonal component in the inertia tensor. This is 
the z component of the term describing the mixing of the rotational and isorotational 
momenta. This term is almost constant. For the hedgehog it has to be there to remove 
the symmetry mode from the Hamiltonian. It appears to saturate at a value of about 
73 S.U. for large separations. 

In Fig. H we show the decoupling parameter D, defined in Eq. (j38| ). It is small, 
showing good decoupling. The discontinuity near Q is 0.1 is due to the fact that there 
the harmonic energy becomes zero, and we have mixing with the real zero modes, 
leading to slightly anomalous behavior. This is not important for the quality of the 
solution outside these points, since the conditions for a valley are local. 

In Fig. |9| we show the potential energy V for movement in the attractive channel. 
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Figure 7: The inertial parameters in the attractive channel, obtained by following the 
lowest non-zero mode from the donut. The lower three panels show the x, y and z 
components of the (iso)rotational moments of inertia (from bottom to top). Here solid 
lines represent rotational moments of inertia, dashed lines isorotational and dashed- 
dotted lines the off-diagonal component of the inertia tensor coupling the two. The 
upper panel shows the translational mass. The solid lines give the xx-, the dashed 
lines the yy- and the dashed-dotted lines the zz-component. Any component of the 



zero-mode inertia-tensor not shown is identically zero. Q is defined in Eq. (34). 
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Figure 8: The decoupling measure D, Eq. fl38|), in the attractive channel, obtained by 
following the lowest mode. Q is defined in Eq. ([34]) . The violent jump is due to poor 
convergence when the RPA frequency of the mode followed goes through zero. 
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Figure 9: The potential energy for movement in the attractive channel, obtained by 
following the lowest mode. Q is defined in Eq. fl53|). 
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Figure 10: The size of the solution in three directions in the attractive channel, obtained 
by following the lowest mode. Q is defined in Eq. fl34l). The solid line gives R, the 
dotted curve is J (x 2 ), the dashed curve J (y 2 ), and the dash-dotted curve J (x 2 ). 



Note that we have defined Q such that the mass is exactly 1 so this parameter does 
not play a role here. 

In Fig. |l(] we show the size of the solution. We also exhibit the coordinate R, 
defined in Eq. ([H]). This can be used to remap all result in terms of this coordinate. 
In particular, one finds that M(R) = (dQ/dR) 2 , all other quantities being scalars. 



5.2 The hedgehog and the hedgehog- hedgehog channel 

Let us first study the hedgehog solutions, to see whether there is a way out of the 
dilemma that the hedgehog is not connected to the hedgehog channel. As expected we 
find two sets of three negative energy modes, one of which corresponds to the hedgehog 
coming apart by an isorotation of the individual hedgehogs, and the other to the fission 
of the hedgehog into two hedgehogs. Since we have three axes, both sets of modes are 
threefold degenerate. In our numerical work following one of the fissioning modes, we 
had problems due to a mode that appeared spurious for the hedgehog (and was thus 
discarded), which becomes non-spurious for finite separation. Similar things always 
occur when we approach a state where an additional symmetry is realized, such ad 



in nuclear physics, if we approach a state of axial symmetry [36]. In these cases the 



RPA frequency has a well defined limit as one approaches symmetry, since B and V 
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Figure 11: A schematic representations of the eigenvectors of the RPA problem around 
the B = 2 hedgehog. The arrows in the box denote the direction of change of the given 
component of the eigenvector. The four rows of three points in the left half of the box 
denote from bottom to top, the x, y and z components, with the upper line giving the 
change in 0j. The two uneven rows on the right denote the change in A (lower) and T 
(upper). We have also indicated the squares of the harmonic frequencies. 
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(short-hand for the covariant second derivative of V) approach zero at the same rate. 
Since the harmonic frequency is the square root of V/B, we cannot have that V is 
finite while B goes to zero. Actually as discussed in great detail in appendix [B[ this 
case is exactly realized in our numerical calculations! We find it hard to believe that 
the Skyrme model itself is so pathological, so we must assume that this is a problem of 
the AM Ansatz. That we do not see a way to smoothly connect the pole configuration 
in the hedgehog-hedgehog channel to the hedgehog is probably closely tied to this 
problem. It is worrisome and probably signifies that for strongly interacting Skyrme 
systems at small distances we should not take the AM Ansatz to be the final word. For 
the NN force this may not be so important, since the Skyrme model itself is suspect 
at these distances. 

Studying the hedgehog channel is no major problem. The technique we used is to 
start with a pole configuration obtained from the "turning-point approach" discussed 
in Sec. ^ and iterate until convergence. At the point thus obtained we follow the 
relevant mode in both directions to obtain a picture of this channel. Fortunately the 
initial guess was quite close to the actual path. 

Some of the results are presented in Figs, p"2|-15. In Fig. 12 we exhibit the inertial 
parameters for rotations, isorotations and translations (i.e., the mass) as a function 
of the coordinate Q. Since the configurations have axial symmetry around the axis 
of separation (x-axis) the yy and zz components of all these tensors are equal. The 
mass is almost constant since the energy is almost constant, and the breaking of Galilei 
invariance remains a minor nuisance. The moment of inertia for rotations increases 
drastically as Q increases, for rotations orthogonal to the axis of separation (x-axis). 
Parallel to this axis it decreases as expected. 

There is quite some structure in the off-diagonal terms. The inertia for ordinary 
rotations around the x-axis is identical to that for iso-rotations around the same axis, 
whereas the mixing term is the opposite of these values. This of course shows that the 
states are invariant under the action of simultaneous rotations and iso-rotations along 
this axis. Quantum-mechanically this implies that there can be no dynamical motion 
along this axis, and thus no kinetic energy proportional to (I x + J x ) 2 . Furthermore 
there is an (identical) off-diagonal term connecting the y and z components of rotations 
and isorotations. 

In Fig. [13] we show the decoupling parameter D, defined in Eq. (|38|). This parameter 
is extremely small everywhere, showing excellent decoupling, i.e., one coordinate is 
enough to describe the dynamics in this channel. The small discontinuity near Q = .1 
is probably due to some solutions that are slightly less well converged than at the 
neighboring points. This is of course of no relevance, since the discontinuity in D is so 
small. 

In Fig. n we show the potential energy V for movement in the hedgehog-hedgehog 
channel. This shows the expected repulsion in this channel. Note that the energy of 
the B = 2 hedgehog is 219.8 S.U., considerably higher than what we have reached in 
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Figure 12: The inertial parameters in the hedgehog-hedgehog channel, obtained by 
following the lowest mode with the relevant symmetries. The lower two panels show 
the x and the y = z components of the (iso) rotational moments of inertia (from 
bottom to top). Here solid lines represent rotational moments of inertia, dashed lines 
isorotational and dashed-dotted lines the off-diagonal component of the inertia t ensor 
coupling the two. The upper panel shows the translational mass. The solid lines give 
the xx-, the dashed lines the yy- and the dashed-dotted lines the ^-component. Any 
component of the zero-mode inertia-tensor not shown is identically zero. 
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Figure 13: The decoupling measure D, Eq. fl38f), in the hedgehog- hedgehog channel, 
obtained by following the lowest mode with the correct symmetry. 
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Figure 14: The potential energy for movement in the hedgehog-hedgehog channel, 
obtained by following the lowest mode. 
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Figure 15: The size of the solution in three directions in the hedgehog-hedgehog chan- 
nel, obtained by following the lowest mode. Q is defined in Eq. ([34]). The solid line 

gives i? 2 , the dashed line J (x 2 ), the dashed-dotted line J (y 2 ) = J (z 2 ). 



the current calculation. 



In Fig. 15 we show the size of the solution. We also exhibit the coordinate R, defined 



in 21. One can see that the transverse size of the solution shrinks as R decreases. 



5.3 The repulsive channel 

A study of the repulsive channel appears to be more challenging. The first problem 
is to find the correct mode to start from. In order to find this solution we vary the 
parameter D, the size of the triangle, Eq. (p0|), for a given value of the parameter 
A2 = <W — ^t- For a suitably large value of A (i.e., large separation) we then look for 
those values of D such that the valley conditions are satisfied. The problem arises that 
there are at least two solutions that look like they could evolve to a product Ansatz 
solution for large separation. We have followed both of them and rejected one which 
appears to evolve to a hedgehog that is too large. A further problem is that in this 
case the relevant mode is always the highest in RPA frequency. Since there is a mode 
the preserves the symmetry lower in the RPA spectrum, this violates the adiabatic 
assumption, and leads to instabilities. The other candidate appears to describe a 
solution that is better behaved. The hedgehogs seem a bit small, but we cannot follow 
to arbitrarily large separation. 
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Figure 16: The inertial parameters in the repulsive channel, obtained by following the 
lowest mode. The lower two panels show the x and the y = z components of the 
(iso)rotational moments of inertia (from bottom to top). Here solid lines represent 
rotational moments of inertia, dashed lines isorotational and dashed-dotted lines the 
off-diagonal component of the inertia tensor coupling the two. The upper panel shows 
the translational mass. The solid lines give the xx-, the dashed lines the yy- and the 
dashed-dotted lines the zz-component. Any component of the zero- mode inertia-tensor 
not shows is identically zero. 
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Figure 17: The decoupling measure D, Eq. (|38D, in the repulsive channel, obtained by 



following the lowest mode with the correct symmetry. 



The inertia parameters, Fig. [TBI, show similar behavior as in the previous case. The 
inertia for ordinary rotations around the x-axis is identical to that for iso-rotations 
around the same axis, whereas the mixing term is the opposite of these values. Like in 
the hedgehog-hedgehog channel this implies the absence of terms of the form (I x + J x ) 2 
in the Hamiltonian. There are no further off-diagonal terms. 



In Fig. 17 we show the decoupling parameter D. The decoupling measure is large, 
especially so at small separations. This indicates that a one-dimensional approach to 
the repulsive channel is not sufficient. There are two modes that contribute significantly 
to D. The manifold of these three modes corresponds exactly to the parametrization 
of Eq. (|2ll|). The one mode that preserves the reflection symmetries is always higher 
than the mode we follow, and the one that breaks the symmetry considerably lower in 
the spectrum (in most cases very close to zero, showing that the energy is almost flat 
in this direction), which has to do with being at a maximum of the energy (cf. Fig. |j). 
We should at least include all three these coordinates to get a reasonable result for a 
decoupled manifold, and thus for the collective Hamiltonian. The oscillations in D are 
due to problems to bring the solutions to convergence. For such a large decoupling 
measure, every minor flaw is amplified significantly. 

In Fig. [18] we show the potential energy V for movement in the repulsive channel. 
The repulsion is indeed even stronger than in the hedgehog-hedgehog channel. 
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Figure 18: The potential energy for movement in the repulsive channel, obtained by 
following the lowest mode. Q is defined in Eq. (El) 




Figure 19: The size of the solution in three directions in the repulsive channel, obtained 
by following the lowest mode. Q is defined in Eq. fl34|). The solid line gives R, the 



dashed line \/(x 2 }, the dashed-dotted line J \y 2 ) = J(z 2 ). 
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Figure 20: The potential energy as a function of R. The solid line gives the result for 
the attractive channel, dashed for the hedgehog-hedgehog channel, and dash-dotted 
for the repulsive channel. 



5.4 NN projection 



As discussed in great detail in Ref. |HJ and Refs. [§, [TIJ, one needs at least the 
three channels discussed here to project out a nucleon-nucleon force from the clas- 
sical Skyrmion-Skyrmion interaction. The standard trick is to use the parametrization 



V(R) = vi(R) + v 2 (R)W + v 3 {R)Z, 



(35) 



where W and Z take the values —1,-2, —1, 4 and 3, in the attractive, repulsive and 
hedgehog-hedgehog channels, respectively. One then associates quantum operators 
with W and Z, which allows for the evaluation of a NN projection of V and by using 
the Born-Oppenheimer approximation, the NN-force fpL0| | . Of course this requires 
the identification of the coordinate R along the individual paths, which is probably 
incorrect. A more correct way would be to find the orthogonal geodesic lines that 
couple the points R on one path to the points R! on another, but this appears to be a 
very challenging task (if not outright impossible). We shall first convert the potentials 
presented previously to functions of R. In Fig. |2(] these functions are shown. Of some 
interest is the size of the potential energy in the hedgehog channel, which is much lower 
than in the product Ansatz. 

The mass in each channel is (4§) 2 . 



In Fig. 21 we give In the attractive and hedgehog 
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Figure 21: The mass as a function of R. The solid line gives the result for the attractive 
channel, dashed for the hedgehog- hedgehog channel, and dash-dotted for the repulsive 
channel. 
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channels this is a smooth function, which fact is of course closely related to the quality 
of decoupling. For large R the mass approaches the reduced mass for two hedgehogs, 
36 S.U.. In the repulsive channel the behavior of the mass is surprising, but there 
decoupling is poor as well. One thus finds that the mass to be used in calculations is 
channel dependent, in a way similar to the potential energy. 



6 Conclusions 

It appears that the application of LACM techniques to the Skyrme model is a promising 
path in the study of the NN force. In this paper we have applied these techniques 
to the Atiyah-Manton ansatz, and obtained several interesting results. An analysis of 
the three standard "paths with reflection symmetries" yielded a wealth of information. 
One of the important aspects is that the results for the repulsive channel may not 
be reliable, since decoupling is poor. This should have an impact on the calculation 
of the nucleon-nucleon force from the Skyrme model, but we have not yet considered 
this aspect in the current paper. On the downside it appears that the AM ansatz has 
some undesirable problems (such as the lack of a consistent description of fluctuations 
around the B = 2 hedgehog). It would therefore be preferable to start from an Ansatz 
that is more fundamental than the add-hoc form of the AM Ansatz. 

One such route would be to apply the techniques of this paper to the Skyrme model 
defined on a grid, where the values of U at the grid points are the relevant dynamical 
variables. This would also allow for the addition of a (very physical) finite pion-mass 
term to the Lagrangian. The problem with such an approach is that, instead of a 
few discrete modes, one obtains a discretized approximation to bands of states (as in a 
solid). Only the band-heads are relevant for the discussion of LACM. The identification 
of the different bands is a complicated and highly non-trivial task. It may actually 
only be feasible for the channels discussed in this paper, where symmetries and their 
representations may be of some help in classifying the different phonon branches. Still 
one would like to be able to connect the different channels, through states without 
the additional symmetry, in order to obtain even more information about the geodesic 
structure of coordinate space (for a pedestrian: about the coordinate dependence of 
the kinetic energy). 

Obviously we have not addressed in any detail the structure of the quantum Hamil- 
tonian, and neither have we used all information derivable from the present calculation 



in the construction of the Hamiltonian. In Ref. |37| a boson model was used to analyze 
the finite N c effective Hamiltonian for the spin and isospin rotational modes of the 
hedgehog. There one concentrated on the leading terms in the large N c limit, which 
happens to be the potential energy. In order to make use of the results of the present 
paper one should extend the work of P?| to include information about the kinetic en- 



ergy. One should then calculate the sub-leading terms in the large- N c limit and match 
those on the ones calculated in the present work. Such an investigation is currently 
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underway. Due to the complexity of this task, this will be described in a separate 
paper. 
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A Reflections and the product Ansatz 

In this appendix we study those solutions of the Skyrme model in the B = 2 sector 
that have a simultaneous reflection symmetry in three orthogonal coordinate planes. 
To that end we assume that it is enough to determine what symmetries the product 
Ansatz has in the limit of large separations, and then require the same symmetry of 
any solution that evolves from it for finite separations. In this limit the U fields of 
the two individual baryon-one hedgehogs commute, and we can study the symmetries 
for a single hedgehog, or the relation of one hedgehog to the other, depending on the 
specific reflection plane. 

The product Ansatz is defined by (we rotate one hedgehog by C and the other by 
C\ to obtain an initial form with maximal symmetry) 

U PA (x\RC) = U^U + = C/(x| - R/2, C)U(x\R/2, C ] ) 

= C[cos/_ + ix- ■ rsin/_]C'C^[cos/ + + ix + • rsin/ + ]C. (36) 

Here we use 

R 

f± = f(x±), x ± = x=f— ei. (37) 

The function / is the chiral profile that describes the hedgehog. 

For future use we introduce the rotation matrix D = D[C] associated with C, 

Dy = ^[Cn&rj]. (38) 

We write 

41 =5 kl {l-25 ki ). (39) 

for the diagonal matrix with one entry — 1, which describes a reflection in a coordinate 
plane, We are looking for those forms of the product Ansatz, where the field U = 
«4 + ir ■ u after a reflection of the coordinate system in any of the three coordinate 
planes is related to itself by a coordinate independent transformation. 

Let us first look at the 2, 3. These transformations do not interchange 

the two hedgehogs, so that we can look for the symmetries of each hedgehog separately. 
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Explicitly we find, using Eq. (^8] 



U + = cos / + + i(Dx + ) ■ t sin / + 

-> cos /+ + z( J Dcr (i) x + ) ■ r sin /+, (40) 

U- = cos /_ + i(D T xJ) ■ t sin /_ 

-> cos/_ + i(£>V i) z_) -rsin/_. (41) 

We require that the components of U± transform as 

ui -> 4, u±-»5 ( V, (42) 

where the matrix is constant. Under the remaining transformation X\ — > —X\ we 
do interchange the two hedgehogs, 

U+ -> cos/_ +i(Dcr (1) L) -rsin/., (43) 
[/_ -> cos/+ + i(I>V 1) x + ) ■ rsin/+. (44) 

Here we require the intertwining symmetry 

u 4 ± -> u* , u± -> S (1) u T - (45) 

The symmetries lead to the conditions 

S®D = D<r®, S®D T = D T a®, £ = 2, 3, 

S®D = D T *W, SWD T = DaM. (46) 

These can converted in equations for the 5"s 

S® = D<r®D T = rF<T®D, i = 2,3, 

S« = D T a^D T = Da^D. (47) 

Now let us see what these equations, which state that we have simultaneous sym- 
metry under all three reflections, have to say about the matrix D. Equating the two 
forms for for £ = 2,3 leads to a condition on D 2 , 

a®D 2 (T® = D 2 . (48) 

The only type of matrix satisfying this condition for % = 2 and 3 simultaneously is a 
diagonal matrix, 

D 2 = diag(±l, ±1, ±1), (49) 

with an even number of minus signs. Thus D 2 is either the identity or a rotation 
around any axis over an angle 7r. The only form for C allowed is then C = exp(z| • n<f>) 
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with = ^ . The second of Eqs. ((47|) shows that = S^ T . We find the following 
three conditions 



sin0n 2 n 3 (l — cos ( 
sin0n^n 2 (l — cos< 
sin (ph\hi (1 — cos< 



0. 

0, 

fii sin ■ 



(50) 
(51) 
(52) 



There clearly exists the trivial solution for sin0 = 0, or0 = kn, with no restriction on 
h. Two more interesting classes of solutions exist as well, where can be any multiple 
of vr/2 



n = (1, 0, 0),n = (0, cosa, sina). 
Let us now calculate the matrices S for the interesting classes of solutions. 

1. = 0: (HH-channel) 

2. h — (1,0,0), = ±7r/2: (repulsive channel) 

S«=diag(-1,-1,-1), 

5( 2 )=diag(l,l,-l), 
5^=diag(l,-l,l). 

3. n = (0, cos a, sin a) (sin0 = 1 only): (attractive channel) 



(53) 



a 



(i) 



S (2) 



S (S) 



( cos 2a 2 cos 2 a sin a 2 cos a sin 2 a 

2 cos 2 a sin a 1 — 2 cos 4 a —2 cos 3 a sin a 
\ 2 cos a sin 2 a —2 cos 3 a sin a cos 4 a + sin 4 a 

-cos 2a —2 cos 2 a sin a —2 cos a sin 2 a 

-2 cos 2 a sin a cos 4 a + sin 4 a —2 cos a sin 3 a 
-2 cos a sin 2 a —2 cos a sin 3 a — 1 + 4 cos 2 a — 2 cos 4 a J 



(54) 

(55) 

(56) 
(57) 

(58) 
(59) 

(60) 



The standard reflection symmetries are usually taken to be the case = and = 7r/2 
with h = (1, 0, 0) and (0, 0, 1). For this last unit vector we have 



5-d) 

5 (2) 



a 



(i) 



diag(-l,l,l), 
diag(l,l,-l). 



(61) 
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B Analysis of the behavior near the B = 2 hedgehog 



Here we shall analyze the behavior of the mass matrix and the potential energy near 
the B = 2 hedgehog. Starting from 

1 A 1 



p(x,t) 



+ 



+ 



(62) 



x 2 + (t-T) 2 x 2 + t 2 x 2 + (t-T) 2 ' 

we shall study the effect of changing the parameters Aj and Tj. The advantage of such 
a choice is that the matrix A 4 is Abelian, and for any change the U field is of hedgehog 
form, 

U = exp(ir ■ x6{x)). (63) 

It is also not a restriction to study only this limited subset of modes, since in a numerical 
calculation we find that they decouple from the remaining modes. 

To calculate the (9-field it is only required to evaluate a simple integral, 

/oo 
f(x,t)dt, 
-oo 

f(x,t) = d x 2 hip = d x2 p/p. (64) 

We shall use the subscript to denote the default choice of JNR parameters of Eq. ( |6"2"D . 
If we now analyze the variations of /, we find that to lowest order 



f = fo + 
5Xo 



1 

Po 



5Ai 



1 



x 2 + (t + T) 2 \x 2 + (t + T) 2 
1 



-fo + 



X 2 +t 2 

SX 3 



x 2 + t 2 



fo) 



+ 



1 



x 2 + (t - T) 2 
2(t + T)8Ti 



x 2 + (t- T) 2 
{ 2 



fo 



{x 2 + (t + T) 2 ) 2 {x 2 + (t + TY 



fo + 



2t5T 2 



2\2 



(x 2 + t 



2(t-T)ST 3 



x 2 + t 2 



fo + 



{x 2 + (t-T) 2 ) 2 \x 2 + (t-T) 2 



-fo 



(65) 



From this we can easily read off that there are four out of a total of six possible changes 
that leave 9 invariant. These are 

1. (5Ai, 5A2, 5A3) = e(l, A, 1). In this case the three first terms inside the curly brack- 
ets combine to edpo, and the three second terms combine to — epo/o = — e^Po- 
These two contributions cancel. Of course this is the simultaneous rescaling of 
all three A's, which is a zero mode for all pole configurations. 
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2. (SXi, SX2, SX3) = e(l, 0, —1). The reason here is time reversal invariance; thus 



/p x 2 + (t + T) 2 {x 2 + (t + T) 2 fo } dt 
is invariant under the interchange T — > — T, and the two contributions cancel. 

3. (ST 1 ,ST 2 ,6T 3 ) = e(0, 1,0). The integrand is odd under the interchange t — > —t, 
and the integral vanishes. 

4. (5T 1; 5T 2 , 5T 3 ) = e(l,0, 1). The time-odd parts of the integrand give 0, and the 
two time-even parts cancel. 

Of course the modes listed as 3 and 4 can be combined to the other trivial zero, a 
shift in the time-origin, which exists for all possible pole configurations. For the case 
of the hedgehog we are thus left with two additional unphysical modes, which are due 
to some residual gauge invariance [17|, which is unbroken only for the particular choice 
of poles discussed here. 

The mass matrix, which can be directly related to the first derivative of U w.r.t. 
the parameters must have zero eigenvalues in this space. It is not so obvious, and we 
shall argue also not true, that V has zero eigenvalues in the same subspace. If this 
were true we would feel happy discarding this part of the space, since no dynamical 
information is contained in these modes. 

To study this problem we look at both B and V in the basis (components, in order, 
8X1, 5X 2 , 5X 3 , 5Ti, ST 2 , 8T 3 ) 

e a = 1 (1/A,1,1/A, 0,0,0), 

yfi + 2/A 2 

e 2 = -^=(0,0,0,1,1,1), 

e 3 = -^=(1,0,-1,0,0,0), 

e 4 = . 1 (1, -2/A, 1,0, 0,0), 

e 5 = -^=(0,0,0,1,0,-1), 

e 6 = -^(0,0,0,1,-2,1). (66) 
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In this basis B and V have the six-by-six matrix form 

/ 





B = 






Vo 
















77.593 
17.385 








17.385 
3.895 




\ 





07 



(67) 



V 



(68) 



/ \ 



1.175 .5474 

24.176 5.426 

5.426 1.212 

V .5474 .2464 J 

If one compares these two matrices, the problem is obvious: the two-by-two block 
spanned by the third and sixth vector is zero for B (as we had argued before), but 
it is non-zero for V. Actually the V block has one zero eigenvalue, but there is one 
eigenvector with infinite harmonic frequency (u 2 = V / B)\ In the numerical calculations 
we see this eigenvalue coming down to finite values as we move away from the hedgehog. 
We believe that this pathology has nothing to do with the Skyrme model, but is 
caused by the Ansatz. The other two-by-two block (fourth and fifth components) 
has an eigenvalue very close to zero, but the same eigenvector and eigenvalue appears 
both in V and B. This reflects the fact that the chiral profile is, for large A and T, 
approximately a function of A/T 2 . This is only troublesome for numerical calculations, 
and causes no theoretical problems. 



C Self consistent solutions in ID LACM 

Suppose we start from a point £° where 



(69) 
(70) 



are simultaneously satisfied. We could also calculate the left eigenvectors g a of the 
RPA matrix, which we assume to be normalized as 

Cf% = (71) 

We now crank by adding g\ multiplied by a small number to £o to obtain a new point 
£i where these equations are not longer satisfied. More specific, 



N 



(72) 
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Va-Mi^i^f^) = («<">) 2 #tfi), (73) 
where we assume that by taking a small step from a point where the equations are 
satisfied, we have Ai 3> A M . 

(This last condition is only schematic. Since A^ form the components of a vector, 
we should require that 

B ll \\ > B^Aj. (74) 

Here we have used the fact that B is diagonal.) 

Using a linear approximation to the first derivative of V, 

V a (^i + AO = V a (^) + V^i&Af, (75) 

and assuming that the RPA eigenvectors do not change (which is not right, but iterative 
use of our update algorithm should be able to correct for this), we find that the choice 

JV 

WGOA^ = - £ (76) 

fj,=2 

will lead to a point where the desired equations are satisfied. (We have neglected one 
set of first order contributions, originating in the change of the eigenvector f 1 . This 
corresponds to assuming that the derivative of the RPA matrix is 0.) Let us now try 
to simplify the solution to (f76|) as much as possible: 

V a -,pB^ = £/>^) 2 ^, (77) 



B aP {V7^ = Ef^ i>M) r 2 9l, (78) 

= E&tfwr 1 ^- ( 79 ) 

Here we have used the fact that in practice V is diagonal (in general we can always 
make a choice of basis such that this happens). If we now use A^ = V a g a , we have 

N 
li=2 fx' 



N 
N 

-Y.9l{V^)- l 9lV rr (80) 

H=2 
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